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Abstract 

We introduce a theory that exposes the fundamental and previously overlooked connection be¬ 
tween the correlation among electrons and the degree of quantum coherence of electronic states in 
matter. For arbitrary states, the effects only decouple when the electronic dynamics induced by 
the nuclear bath is pure-dephasing in nature such that [H ^= 0, where is the electronic 
Hamiltonian and is the electron-nuclear coupling. We quantitatively illustrate this connection 
via exact simulations of a Hubbard-Holstein molecule using the Hierarchical Equations of Motion 
that show that increasing the degree of electronic interactions can enhance or suppress the rate of 
electronic coherence loss. 
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Understanding the behavior of electrons in matter is fundamental to our ability to char¬ 
acterize, design and control the properties of molecules and materials BE]. Electronic 
correlations M and decoherence I3EI are two basic properties that are ubiquitously used 
to characterize the nature and quality of electronic quantum states. Correlations among 
electrons arise due to their pairwise Coulombic interactions, that lead to a dependency 
of the motion of an electron with that of other surrounding electrons. These correlations 
determine the energetic properties of electrons in matter and the character of their energy 
eigenstates mu- In turn, decoherence in molecules typically arises due to the interactions of 
the electrons with the nuclear degrees of freedom [TUHT^ . The nuclei act as an environment 
that induces a loss of phase relationship between quantum electronic states. Establishing 
mechanisms for electronic decoherence is central to our understanding of the excited state 
dynamics of molecules [T3HT6], to the development of useful approximations to model corre¬ 
lated electron-nuclear dynamics [niiiH], and to the design of strategies to preserve electronic 
coherence that can subsequently be exploited in quantum technologies mm- 

While electronic correlation and decoherence have been amply investigated separately, 
the connection between the two, if any, is not understood. This is partially due to the fact 
that usual dehnitions of electronic correlation, such as correlation energy [21] or natural 
occupation numbers [22] , are only applicable to pure electronic systems [231123] and do not 
allow addressing this fundamental question. For this reason, it is unclear if decoherence 
can induce changes in correlation and, conversely, if correlations can modify the coherence 
content of a quantum state. 

Here we demonstrate that electronic correlation and decoherence are coupled physical 
phenomena that need to be considered concurrently. We do so by extending the concept 
of electronic correlation to open non-equilibrium quantum systems, and showing that elec¬ 
tronic correlation modulates the degree of entanglement between electrons and nuclei, and 
thus the degree of electronic decoherence. Conversely, we also show that the electronic de¬ 
coherence modulates the degree of electronic correlation, as evidenced by the correlation 
energy. Further, we isolate conditions under which electronic correlations and decoherence 
can be considered as uncoupled physical phenomena and show that they are generally vi¬ 
olated by molecules and materials, demonstrating that the connection between electronic 
correlation and decoherence is ubiquitous in matter. These formal developments are quanti¬ 
tatively illustrated via numerically exact computations in a Hubbard-Holstein molecule that 
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show that increasing the electronic interactions can strongly modulate the rate of electronic 
coherence loss. 

To proceed, consider a pure electron-nuclear system with Hamiltonian "H = + 

where is the electronic Hamiltonian, the nuclear component, and the 
electron-nuclear couplings. Here, is dehned as the residual electron-nuclear interactions 
that arise when the nuclear geometry deviates from a given reference conhguration (e.g., the 
optimal geometry). The electronic Hamiltonian = Hq + can be further decomposed 
into single-particle contributions Hq (e.g. Hartree-Fock) and residual two-body terms . 
The latter arise from Coulombic interactions that cannot be mapped into one-body terms 
and introduce correlations among the electrons. The associated non-interacting Hamiltonian 
is obtained when = 0, and is given by "Ho = Hq + . 

To extend the concept of electronic correlations to open non-equilibrium quantum sys¬ 
tems, we require a correlation metric and a reference uncorrelated state for each electron- 
nuclear state. To construct the reference state, we imagine a hctitious process where for 
each physical time t the term in the Hamiltonian is turned off adiabatically slow along 
a hctitious time coordinate r (see Fig. SI in the Supporting Information (SI)). Specihcally, 
we suppose that the Hamiltonian of the system is of the form 

(e>0), (1) 

where the second term is considered as a perturbation to the "H-induced evolution. The 
physical evolution along t occurs at the r = tq —)■ —oo limit of the (t, r) space for which the 
Hamiltonian is in its fully interacting form HeiTo) = H. In this limit, the state of the fully 
interacting system is given by 


pit) = '^(^iit)a*it)\'ipi){'ipj\, ( 2 ) 

where 1'^*) are eigenstates of H i'H\'ipi) = Ei\'ipi)). The uncorrelated reference state is gener¬ 
ated by adiabatically turning off, in the Interaction picture, the term in the Hamiltonian 
in the r = tq to r = 0 interval, i.e. 


= liin lim f/,/(0,ro)p(t)Hj/0 ,Tq), (3) 

e^O TQ^—oo 

where t') is the evolution operator in Interaction picture. The latter is dehned by the 
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Dyson series [1] f/e/(r, r') = I + where 

j dTne~^^^^^Vl{Tn)U^j~^\Tn,T'), 

Vj{t) = —Uq{t)V^Uo{t) is the —V^ operator in Interaction picture, and Uo{t) = is 

the perturbation-free evolution operator. 

Equation (|^ captures changes in p{t) that are generated by the process of turning off 
in the presence of a nuclear environment. It has the desirable property that p'^{t) = p{t) when 
= 0, and it reduces to the usual adiabatic connection for isolated electronic systems when 
= 0. Note that we have chosen instead of the full evolution operator U{t) = 

Uq{t)U^i{t) to generate the uncorrelated states. This is because the Uq{t) component of 
U{t) leads to changes in p{t) due to electron-nuclear entanglements that are present even 
when = 0. By contrast, U^i{t) solely captures electron-nuclear entanglements that can 
be modulated by the electron-electron interactions. 

Switching off interactions adiabatically generates exact eigenstates of the non-interacting 
system from those of the interacting system via the Gell-Mann and Low theorem (GMLT) [U 
m. The GMLT states that given an eigenstate l-^j) of the interacting T-L, if the limit 

liin 10^) = -oo)|V’i), (4) 

(where Ai = (' 0 j|f/e 7 (O, —cxD)|' 0 j)/|(' 0 j|f/e/(O, —cxD)|' 0 j)| and = 1 because the | 0 j) are 

chosen to be normalized) exists, then lime_,.o \4>j) = \4>i) is an eigenstate of the non-interacting 
■Ho- Applying the GMLT in Eq. (|^ we arrive at the uncorrelated reference state that 
corresponds to p{t) in Eq. (|^, 

p"(t) = (5) 

ij 

Here, we have assumed that lime_>o AiA* = exists even when the phase factors 

Ai ~ ce are known to be ill-behaved as e —)■ 0 [1]. While the A^ introduce convergence 
issues at the wavefunction level, observable quantities, including the density operator, should 
remain hnite during the unitary evolution. 

As a physical measure of electronic correlation in electron-nuclear systems we choose the 
energetic difference between the correlated and uncorrelated state: 

Ecor(t) = Tr[p(f)H] - Tr[/)"(f)Ho]. ( 6 ) 
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This quantity measures energetic changes in the electron-nuclear system that are introduced 
by the process of turning off during the adiabatic connection in Eq. ([^, and parallels 
a common metric for correlation [2Tj used in closed electronic systems. Note that any 
energetic measure of correlation based on the properties of the electronic subsystem alone 
is not appropriate since it will unavoidably include relaxation channels due to interactions 
with the bath. Further note that definitions of correlation based on the non-idempotency 
of the single-particle electronic density matrix |22l |25] are not applicable since the non- 
idempotency can arise due to correlation or due to decoherence [2S]- (see Ref. [21112H] for 
measures claimed to operate in open quantum systems). 

As a basis-independent measure of decoherence we employ the purity P{t) = Tr[/3g(t)] 
where Pe{t) = TTB[p{t)] is the A^-body electronic density matrix obtained by performing a 
partial trace over the nuclear bath. The purity P = 1 for pure states and P < 1 for mixed 
states. For pure electron-nuclear systems, the decoherence of the electronic (or nuclear) 
subsystem is solely due to electron-nuclear entanglement. Thus, in this regime, the decay of 
P also measures the degree of electron-nuclear entanglement. 

In this context, it is now readily seen why correlation and decoherence are strongly 
connected. For this, first note that the coherence content of pe(t) and p^it) are gener¬ 
ally different. To see this, consider p{t) = |T(f))(T(f)| in Eq. @ for which |T(f)) = 
Q;j(f)|'0j). In light of the Schmidt decomposition [20], |T(t)) can be written as |T(f)) = 
\/AWki(f))l^i(^)), where |sj(f)) and \bi{t)) are, respectively, orthonormal electron and 
nuclear states, and a/\ are the Schmidt coefficients (Xli= 1, A, > 0). In the Schmidt 
basis. 


p{t) = ^/xY\si)\bi){bj\{sj\. (7) 

In terms of {Aj}, the purity of the electronic (or nuclear) subsystem is P{t) = 

In turn, the uncorrelated state p'^{t) = |<h(t))(<I)(f)| (Eq. 0 ) is associated with |<h(f)) = 
Q;j(f)e*®*|0j). Under the Schmidt decomposition, |<h(t)) = ^/pi{t)\Si{t))\Bi{t)) and 

the resulting purity is P“(t) = Since |‘h(t)) 7^ |^(^)), fhe set {pi\ is different from 

the set {Aj} and therefore the purity for the correlated state and its reference uncorrelated 
counterpart generally differ. That is, for 7^ 0, modulates the degree of coherence of 
electronic states. 
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Consider now the influence of on the correlation energy [Eq. (|^], 

Ecor{t) = ^ \ai{t)\‘^{Ei - Si), 
i 

where = Ei\'ipi) and = Si\(j)i). For ^ 0, E^or will change if El^^ changes 

because Ei and Si vary differently as is modified. That is, influences Ecor because 
it modulates the response of the electron-nuclear system to . 

Decoherence and correlation decouple when 

= 0, (9) 

for 7 ^ 0. When Eq. (|^ holds, the does not introduce electronic relaxation and the 
system-bath dynamics is pure dephasing. To see how this sufficient condition arises, consider 
the decoherence case first. For the purity of pe{t) and p'^{t) to coincide, the evolution operator 
in Eq. (|^ must not change the degree of entanglement between electrons and nuclei. For 
this to happen, f/^/ must be of the form 

f/e/(0,-cx)) = f/5(0,-cx)) 0 I^, (10) 


where is a purely electronic operator and is the identity operator in the nuclear Hilbert 
space. Under these conditions, and in light of Eq. 

where |s')(s' | = lime^oU5(0, -oo)\si){sj\U^j\0, —oo). Since the Schmidt coefficients for 
p'^{t) are the same as those of p{t) (c/. Eq. ([^) the purity of the two states is identical. 
For U^i to be of the form in Eq. (10), Vi{t) must be a purely electronic operator, i.e. 
Vj{t) = 0‘^(r)®I'®, where is an operator in the Hilbert space of the electronic subsystem. 
This is guaranteed when Eq. ([^ is satisfied. Specifically, 

where we have used the fact that = 0 and the condition in Eq. ([^. We arrive at 

the desired form 

Vi{t) = = ( 11 ) 

by taking into account that [U'^, TT'®] = 0, and the fact that [U'^, = 0 for Coulombic 

systems since and are both functions of the position operators. 

The correlation energy also becomes independent of when the commutation relations 
in Eq. (9) are satisfied. To show this, we contrast Ecor with the correlation energy Ecor that 
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would have been obtained if is not allowed to influence the response of the system as 
is adiabatically turned off in Eq. ([^. Specihcally, 

£«>(() = H\pm] - (12) 

where the reference state = hmT-o^._oo,e^.o \u'^j{0,To)p{t)Ul\{0,TQ)\ is obtained by 

setting =0 throughout the adiabatic process, i.e. —oo) = U^i^O,— oo)\hsb^q. 

The interaction potential V/(r) in —oo) is given by 

V}'(r) = Vi{t)\hsb^o = ® (13) 

where we have used the fact that [H^,H^] = [V^,H^] = 0. If Ecor{t) = Ecol{t) the 
correlation energy is independent of For this to happen, the identity Vj (r) = V/(r) 

must be satisfied such that and p'^(t) coincide. Since V/(r) is identical to the limiting 

V7(r) in Eq. (11), by the same argument employed to arrive at Eq. (0 we conclude that 
Ecoiit) = Ecol{t) when Eq. (9) is true. 

From the perspective of the correlation energy, when Eq. (|^ is satisfied Ecor is purely 
determined by the electronic subsystem. This is because + H^) remains constant as 

is turned off adiabatically (as can be seen by writing the Heisenberg equations of motion 
for + H^). From the perspective of the purity, Eq. ([^ guarantees that the effect of 
the bath will be the same for the correlated system and its uncorrelated counterpart, thus 
eliminating a possible dependence in the decoherence dynamics. Note that even for 
stationary Born-Oppenheimer (BO) states it is not possible for decoherence and correlation 
to be uncoupled unless Eq. (|^ is satished. This is because even when stationary BO states 
are not entangled, the corresponding uncorrelated state generally will be. 

The pure dephasing condition [Eq. (|^] is generally violated by molecules and materials, 
indicating that the connection between electronic correlation and decoherence is ubiquitous 
in matter. Nevertheless, pure dephasing dynamics can arise when the frequencies associated 
with nuclear motion are far detuned from the electronic transitions such that the nuclear 
dynamics does not lead to electronic transitions in the correlated and uncorrelated system, 
as can be the case in semiconducting quantum dots [221 ED]- Under such conditions, ^ 
® l-Un) (U„|, wliere {\En)} are the eigenstates of , and the En are nuclear operators 
dehned such that [V^, = 0. 

We now quantitatively illustrate this connection using a neutral two-site, two-electron, 
Hubbard-Holstein model with zero net spin as an example |1]; a minimal molecular model 
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that violates the commutation relations in Eq. (i and satishes = 0 as is expected 

for molecules. Here the electrons are described by the Hubbard Hamiltonian 


= -to ^ {d\^d2a + dladia) + U + U2th2^) (14) 


where d]^ (or dia) creates (or annihilates) an electron on site i with spin a and satishes the 
usual anticommutation relations {di^, dj^,} = The quantity hjo- = dj^di^ is the num¬ 

ber operator, to is the hopping parameter, and U is the energy penalty for having two elec¬ 
trons on the same site. The Hubbard Hamiltonian can be decomposed into a Hartree-Fock 
component = -toJ2ae{^,i}idiad2^ + dladia) + 2 t/_<^) - U Y,i,A^io){ni-a), 
and a two-body term — Hq , where the expectation value (hjo-) = 1/2 is over the 

equilibrium thermal state. The nuclei are described as four baths of iV™ harmonic oscillators, 
with Hamiltonian 




Nr 


EE 

m=l j=l 



1 2 

+ 



(15) 


where Xmj is the mass-weighted displacement away from equilibrium for the jth harmonic 
oscillator in the mth harmonic bath, p^j is the momentum conjugate to Xmj and Umj its 
oscillation frequency. We assume that each set of harmonic oscillators couples to an inde¬ 
pendent electronic conhguration of zero net spin. Specihcally, we choose 


= Fihith^ -f F2h2tn2i + ^ 3711 ^^ 2 ^ + F4n^h2^, (16) 

_. _/V^ 

where F^ = X]j=i ^mjXmj is a collective bath coordinate of bath m. The effective electron- 
nuclear coupling is specihed by the spectral density Jm{^) = f ~ ^mj) which 

is assumed to be same for all the states and of Debye form J{(jj) = t] ^ 2 ■ Here 7 is the 
characteristic frequency of the bath and the parameter p effectively determines the electron- 
nuclear coupling strength. 

The electronic dynamics generated by this model is propagated exactly using the Hi¬ 
erarchical Equations of Motion approach [HTHM], a non-perturbative and non-Markovian 
theory of reduced system dynamics. As an initial state, we consider a separable electron- 
nuclear state p{0) = Pe(0) ® Pn(0) where the nuclei are initially at thermal equilibrium 
Pn(0) = exp(— /9if^)/TrB{exp(— with inverse temperature (3, and the electrons 
Pe( 0 ) = |4')(\['| in a superposition IT) = ;^(|Ei) -|- IE 2 )) between the ground and hrst 
excited state. 




0 100 200 300 400 500 600 700 800 

W^lh 



20 40 60 80 100 120 140 




FIG. 1: Purity and electronic energy during the evolution of the Hubbard-Holstein 
molecule (/d = 1/to, hy = O.Sto)- The inset shows characteristic timescales in P{t) obtained 
from an exponential £t P — Pthermai = 6xp(—t/rj), see SI. The dot size measures 

the magnitude of la/ (blue, a, > 0; red, a, < 0; as is small and not shown). Note how 
increasing U can enhance or suppress the decoherence. 


Figure [T] shows the dynamics of the purity and the electronic energy for different elec¬ 
tronic interactions U and effective electron-nuclear couplings y (inset: characteristic decay 
timescales r* in P{t)). The fact that Eq. (|^ is violated is reflected by the energetic relaxation 
of the electrons. The purity observes a sharp initial decay on a ti timescale, followed by a 
slower dynamics on a T 2 timescale that asymptotically leads the electronic subsystem to a 
state of thermal equilibrium. In the presence of electronic correlations, varying U strongly 
modulates the decoherence and relaxation dynamics. By contrast, in the Hartree-Fock ap¬ 
proximation the purity for this model is independent of U and equal to the one for U = 0. 
For U < 3to the decoherence is determined by T 2 - In turn, for U > 4fo the importance 
of T 2 in the dynamics (as characterized by the dot sizes in Fig. is diminished, and the 
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decoherence time is determined by ri. Note how increasing U can enhance or suppress the 
rate of electronic decoherence. Specifically, for rj = O.ltoi increasing U leads to a decrease in 
the decoherence time. By contrast, for rj = 2.0to, increasing U leads to a decrease followed 
by an increase in the decoherence time. As expected, the rate of decoherence is faster in the 
stronger rj case. 

The molecular mechanisms at play in Fig. [T]can be identified by examining the effect of 
changing p and U on the potential energy surfaces (PESs). As detailed in the SI, increasing 
U brings the ground and first excited state closer together in energy, and reduces the differ¬ 
ence in curvature between their PESs. The first effect increases the decoherence rate because 
it increases the nonadiabatic couplings between the two states. Excitation by an incoherent 
bath leads to decoherence [3S] . Thus, the enhanced excitation of the electrons by the thermal 
nuclei increases the decoherence rate. The second effect, by contrast, slows down the deco¬ 
herence. To see this, recall that for a general vibronic state |T) = \Xn) the electronic 

density matrix is given by pe = Ylinrn(Xrn\Xn)\En){Em\- The coherences between states \En) 
and \Em) are thus determined by the nuclear wavepacket overlap Smn = {Xm\Xn) [I2l 1^ . 
By making the PESs look more alike, increasing U slows down the decay of Smn for each 
member of the initial ensemble due to wavepacket evolution in alternative PESs. It is the 
non-trivial competition between these two effects what leads to the intricate dynamics in 

Fig- [3 

Note that the nonadiabatic couplings between the ground (singlet) and first excited 
(triplet) state that are responsible for the first decoherence mechanism arise due to the 
Fshi-f 77 . 24 . F 4 77. 1477.24 terms in . By contrast, the second decoherence mechanism is 

determined by all four terms in and survives even in the absence of singlet-triplet 

couplings. In this limit, increasing U protects the electrons from the decoherence. 

Does decoherence help us reduce the complexity of the many-body electron problem? 
Figure 1^ shows the evolution of the two-particle cumulant (Tr[A 2 ] = Tr[F)F^ — F)F], where 
F)F is the 1-body electronic density matrix) which measures the importance of 2-body 
contributions to pe that cannot be decomposed in terms of T)f [25]. For an uncorrelated 
closed electronic system P^F^ = P^F and Tr[A 2 ] = 0. As shown, instead of reducing the 
complexity, in this case increasing p (and U) enhances the importance of higher order r- 
body electronic density matrices to the BBGKY hierarchy |37j . 

In conclusion, we have shown that the correlation among electrons and the degree of 
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FIG. 2: 


Two-particle cumulant during the evolution in 



quantum coherence of electronic states are strongly coupled in matter. For arbitrary states, 
only when the system-bath dynamics is pure dephasing such that Eq. is satished can 
correlation and decoherence be considered as uncoupled physical phenomena. Investigating 
the consequences of this fundamental, ubiquitous, and previously overlooked connection 
constitutes an emerging challenge in electronic structure and molecular dynamics. 
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The Supporting Information includes plots for the evolution of the electronic density 
matrix, the PESs for representative examples, and a discussion of the mechanisms at play 
in Fig. 
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Understanding the Fundamental Connection Between 
Electronic Correlation and Decoherence: Supporting 

Information 


I. SCHEME OF THE ADIABATIC CONNECTION IN EQS. (l)-(5) 


^ ^ Wo 


GMLT 


, W<(ro) = W 


Arbitrary state p{t’) 


FIG. SI: At each instant of time t, the uncorrelated counterpart p^{t) of a general 
electron-nuclear state p{t) can be obtained by turning off the residual electron-electron 
interactions adiabatically slow via evolution in the Interaction picture along a hetitious 
time coordinate r. 


H. DECOHERENCE DYNAMICS OF THE HUBBARD-HOLSTEIN MODEL 

To clarify the molecular mechanisms at play in the decoherence dynamics of the Hubbard- 
Holstein model shown in Fig. 1, below we discuss the effect of changing U and rj on the 
electronic potential energy surfaces (PESs) and on the dynamics of the electronic density 
matrix. 


A. Dynamics of the electronic density matrix pe{t) 


Figure S2 shows the dynamics of Peit) in the eigenbasis of for four representative cases 
{U = Oto,6to; V = O.lto, 2.0to)- The characteristic relaxation timescale of each matrix 


element was obtained via an exponential £t and tabulated in Table SI As can be seen from 


Figure S2 and Table SI, the most important diagonal elements of pe(t) {pl^(t) and pf'it)) 
decay faster for U = 6fo than U = Oto indicating that increasing the electron-electron 


SI 











interactions generally leads to faster relaxation. For rj = O.lto, the decay of the initial 
coherence between the ground and hrst excited state is largely unaffected by varying 

U. By contrast, for rj = 2.0fo the decay of is slower for U = 6fo than U = Ofo signaling 
that increasing U protects the coherences between these two electronic energy eigenstates. 
Naturally, the thermalization of pe{t) is faster for the stronger electron-nuclear coupling 
ij] = 2.0to) than for the weak electron-nuclear coupling {j] = O.lto)- As described in the 
sections below, these features of the dynamics of pe{t) can be understood by investigating 
the effect of changing U and rj on the PESs, 


Note that the characteristic decay timescales for Pe{t) in Table SI are related to the deco¬ 
herence timescales obtained from the purity (Table [S^. This is because, the thermalization 
of the purity 


F(«) = Tr|p^(i)l = 5^ 


(SI) 




is determined by the individual contributions to the density matrix p^J{t). 


B. Effect of changing U and p on the PESs 


Consider now the effect of changing U and rj on the PESs of the Hubbard-Holstein model. 
For dehnitiveness, we focus on the limiting case where the molecule is coupled to just one 
of the harmonic oscillators in each independent bath. In this case, the Hubbard-Holstein 
model Eqs. (14)-(16) reduces to: 




m=l 


Ti = — to yX (^10-^20- + d^fjdicr) + U+ n2^n2i) + A 


(S2) 


+ CiXifLi^flii + C 2 X 2 n 2 ^h 2 l + -t- CiXihiih 2 ^, 


where the frequency is taken to be at the peak of the spectral density {hw = 0.3to)- The 
strength of the electron-nuclear coupling is determined by {cm}, chosen to be 
While this system does not correspond exactly to the system that is modeled through the 
HEOM approach, it does allow extracting qualitative understanding of the effect of changing 


U and p on the PESs. Figure shows one-dimensional projections of the four PESs for 
this model along X 2 (with xi = x^ = x^ = 0) and xg (with xi = X 2 = X 4 = 0) for different 
choices of U and p. Note that for this simplihed model the PESs along xg and X 4 (or xi 
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FIG. S2: Evolution of the electronic density matrix Pe(t) under the conditions in Fig. 1 with 77 = O.lto 

(top) and r] = 2.0to (bottom). Each panel corresponds to a different matrix element of 

plP{t) = {Ei\pe{t)\Ej) in the eigenbasis of {\Ei)} ordered with increasing energy. Initially 

Pe(0) = |'I')('I'|, where |4') = + 1^2))- Timescales for the decay of pl^it) based on an exponential ht 

are shown in Table [ST| Note that for p = 2.0to the initial coherence between ground and first excited state 

{pl‘^{t)) is protected by increasing U. By contrast, for p = O.lto the decay of pl^it) is independent of U. 
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(a) r]/to = 0.1 
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P 
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pf 
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ff 
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63.45 

-0.0018 

1533 

pf 
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319.5 

-0.0009 
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(b) rj/to = 2.0 
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P 
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r{k/K) 


-0.2633 

6.134 
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3.025 

pf 

0.5289 

1.244 

0.4515 
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pf 
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ff 

0.0372 

2.810 

-0.0166 
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pf 

0.3265 

3.723 

0.1146 

2.616 

pf 
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pf 
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0.0112 

13.75 

pf 

0.0009 

123.8 

0.0033 

23.17 

pf 

0.0565 

2.601 
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pf 

0.0343 

7.220 

0.0023 
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TABLE SI: Characteristic timescales in the thermalization of the density matrix of the 


Hubbard-Holstein model shown in Fig. ^ The data corresponds to an exponential £t of 
the decay — p^J{oo) = pe~^, where r is the timescale and p its weight. 


and X 2 ) coincide. As can be seen in Fig. S3, for weak electron-nuclear couplings {rj = O.lfo) 


the 4 PESs are very similar, while for the stronger rj = 2.0fo the minimum and curvature of 
the PESs generally differ. In both cases, the effect of increasing U is to bring the ground 
and hrst excited state (or second and third excited states) closer together in energy, and to 
reduce the difference in curvature of the PESs associated with the ground and hrst excited 
state. 

The reduction in curvature difference can be quantihed through 




(AF) = (Fi - F2) = 


E« 


-/3en 


(S3) 


which measures the average of the difference in the curvatures Fn = dE^'^\xi)/dxi between 
ground (n = 1) and hrst (n = 2) excited state along a particular nuclear coordinate Xi, where 
(xi) is the adiabatic PES of the n-th electronic state along x*. The average is taken over 
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FIG. S3: Projections of the PESs for the simplihed Hubbard-Holstein model dehned in 


Eq. (S2) along X 2 (with xi = X 3 = = 0), E^'^\x 2 ), and X 3 (with xi = 0:2 = = 0), 

E^'^^Xs), with n = 1, • • ■ ,4. The PESs along X 4 (or xi) are identical to those along 0:3 (or 
X 2 )- Note that with an increase in U the ground and first excited (or second and third 
excited) states come closer in energy. 


the initial nuclear thermal state, where e„ and \xn) are the eigenvalues and eigenfunctions 
of the n-th harmonic oscillator level, and the inverse temperature. As shown in Fig. [S4| 
increasing U generally leads to a decrease in the difference in curvature between the ground 
and first excited state along all nuclear directions. Note that the (AF) for rj = O.lto are ~ 5 
times smaller than those for r] = 2 . 0 fo- 


The decrease in energy difference between the ground and hrst excited state with increas¬ 
ing U, causes the nonadiabatic couplings (NACs), di 2 , between these two states to increase 


(see Fig. S5). The NACs between electronic eigenstates \(j)i{xm)) and \(j) 2 ixm)) are dehned 
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by m 


d 

dl2 = \ {.(l>l{Xra)\^ -|<^2(a:m))| 


(S4) 


Here \4>i{xm)) refers to the i-th Born-Oppenheimer (BO) electronic eigenstate of the Hamil¬ 


tonian (obtained by diagonalizing everything in Eq. (S2) except the nuclear kinetic energy) 


The di 2 measure the coupling between two electronic levels via nuclear motion. An increase 
in the NACs leads to increased excitation of the electrons via nuclear dynamics. As shown 


in Fig. the NACs are ~ 5 times larger for the case of stronger electron-nuclear couplings. 
Further, the NACs increase signihcantly with an increase in U for both values of rj considered 
as a result of the energy levels coming closer together. 


rjlto^OA q!to=2.0 

<AF>(to^'^irt) <AF>(fo^'^/n) 




FIG. S4: Average of the difference in forces (AF) (Eq. (S3)) between the 
excited state for different U and rj along two different nuclear coordinates 
xi = X 3 = X 4 = 0) and X 3 (with xi = X 2 = X 4 ^ = 0). The decrease in (AF) 
in U can increase the lifetime of coherences between these two states. 


ground and hrst 
X 2 (with 
with an increase 


C. Two competing decoherence mechanisms 

As detailed above, increasing U brings the ground and hrst excited state closer together 
in energy, and reduces the difference in curvature between their PESs. As we now discuss, 
these two effects on the PESs lead to competing decoherence mechanisms that underlie the 
dynamics in Fig. 1. 
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FIG. S5: Nonadiabatic 


couplings (NACs) [Eq. (S4)] between the ground (|0i)) and first 


excited BO electronic state (| 02 )) as a function of the nuclear coordinate X 3 (or X 4 ) for 


rj = O.lfo (left) and r] = 2.0to (right). The NACs vanish along the nuclear coordinate xi 


and X 2 - As U increases, |0i) and | 02 ) come closer in energy (see Fig. S3) causing an 


increase in the NACs. 


1. Increasing U decreases the rate of decoherence because it reduces the difference in curvature 
between the PESs 

In pure electron-nuclear systems, decoherence arises because of nuclear evolution in al¬ 
ternative PESs. To see this, consider the electronic density matrix associated with a general 
entangled vibronic state |C(t)) = Yhn \^n)\Xnif))) 

Pe{t) = Trs{|C)(C|} = ^{Xm{t)\Xn{t))\En){Ej\, (S5) 

nra 

where the trace is over the environmental degrees of freedom, the {|E„)} are the eigenstates 
of and the \Xn{t)) is the nuclear wavepacket associated with the n-th electronic state. 
Note that the coherences between electronic eigenstates (the off-diagonal elements in Pe{t)) 
are determined by the nuclear overlaps Snmif) = (Xm(^)IXn(t))- Thus, the loss of coher¬ 
ences in pe{t) can be interpreted as the result of the decay of the Snm during the coupled 
electron-nuclear evolution [T21 ES]- Anything that leads to a decay in the nuclear overlaps 
(anharmonicities in the PES, nuclear motion in high-dimensional space, etc.) leads to deco¬ 
herence. Standard measures of decoherence capture precisely this. For example, the purity, 
the measure of decoherence that we focus on here, is given by 

P{t) = 5^ I (Xm(t) |Xn(^)) P (S 6 ) 
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and decays with the overlaps between the environmental states Snm- 


The effect of increasing U is to rednce the difference in the cnrvatnre between the gronnd 


and first excited state as revealed by (AF) in Fig. S4 This effect is particnlarly important for 
r] = 2.0to while for r] = O.lto the fonr PESs are very similar to one another even for U = Oto- 
Given the initial coherence between the gronnd and first excited state, this redaction in the 
cnrvatnre is expected to slow down the decoherence becanse it leads to a slower decay in the 
overlap of the nnclear wavepackets 512 associated with these two states for each member of 
the initial thermal ensemble. For r] = 2.0to, this feature is clearly reflected in the dynamics 
of Pe^(t) that shows an increase in coherence time from 1.2 to 7.7 to/h as U changes from 


Ofo to 6to (Fig. S2). By contrast, for rj = O.lto, the shape of the PESs is mostly unaltered 


by varying U, hence the reason why in Fig. S2 the pg{t) decay at approximately the same 
rate for U = Oto and U = 6to- 


2. Increasing U increases the rate of decoherence because it reduces the energy difference be¬ 
tween electronic states 


By reducing the energy difference between the ground and hrst excited state, increasing 
U introduces an additional decoherence mechanism in the dynamics that arises because 
the nuclei are initially prepared in a thermal incoherent state. Specihcally, as the energy 
difference between levels is reduced with increasing U, the NACs between such levels increase 


(see Fig. S5). This increase in the coupling leads to an enhanced excitation of the electronic 


degrees of freedom by the nuclear dynamics. Now, excitation of a coherent system by 
an incoherent bath leads to decoherence [35]. Therefore, the enhanced excitation of the 
electronic subsystem by the thermal incoherent nuclear state leads to an increased rate 


of decoherence. For small U, the PESs in Fig. S3 are well separated in energy and this 
mechanism is suppressed. As U increases this mechanism becomes increasingly important 
leading to faster decoherence. 

Note that for p = O.lto this mechanism is expected to be dominant since the PESs in this 
case are essentially parallel. This explains the comparatively long decoherence time observed 
when U = Oto- By contrast, for r] = 2.0to both decoherence mechanisms are expected to play 
a role. This explains why the decoherence is signihcantly faster for rj = 2.0to with respect 
to 7] = O.lto for all U considered. 
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III. TIMESCALES IN THE PURITY DYNAMICS 


U/to=1.5, Ti/to=0.1 U/to=2.0, T|/to=2.0 




tyti lyti 

FIG. S6: Pictorial representation of the timescales associated with the pnrity dynamics for 
r] = O.lfo (left) and r] = 2.0fo (right) is exemplihed in a particular case. For rj = O.lfo, the 
three timescales Ti,r 2 ,T 3 can be associated with the initial decay, growth and final 
thermalization to thermal purity. For r] = 2.0to, Ti captures the initial decay followed by a 
growth/decay captured by timescale T 2 . 


Here we illustrate the meaning of the characteristic decoherence timescales Ti and T 2 in 


Fig. 1 through a particular example. Figure [S6] shows the three timescales ri,r 2 and Ts 
associated with the purity dynamics obtained through a tri-exponential fit. The initial 


decay is captured by ti which has a signihcant contribution to the purity (see Table S2), 
followed by a growth in purity captured by T 2 and subsequently a hnal decay captured by 
T 3 . The contribution of the T 3 timescale is negligible in the r] = 2.0to dynamics and quite 


small for r] = O.lfo (see Table S2), and is not included in Fig. 1. 
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TABLE S2: Characteristic timescales r* in the dynamics of the purity in Fig. 1 extracted 

3 

from a triexponential £t P(t) — Pthermai = exp(—t/rj), where a* quantihes the weight 

i=l 

of each contribution. Note that 03 = 0 for r] = 2.0to and small for 77 = O.lto- 

(a) rj/to = 0.1 (b) rj/to = 2.0 



-^thermal 

2 

y^aiexp{-t/Ti) 

i=l 

u 

ai 

Ti{to/h) 

a2 

T2{to/h) 

O.Oto 

0.846 

0.63 

-0.225 

6.54 

0.5to 

0.830 

0.67 

-0.209 

6.10 

l.Oto 

0.818 

0.72 

- 0.200 

5.71 

1.5to 

0.810 

0.78 

-0.197 

5.25 

2 . 0 to 

0.804 

0.85 

-0.197 

4.73 

2.5to 

0.805 

0.95 

-0.205 

4.09 

3.0to 

0.882 

1.14 

-0.294 

3.00 

3.5to 

0.852 

1.03 

-0.179 

3.07 

4.0to 

0.776 

0.90 

-0.025 

7.67 

4.5to 

0.591 

1.25 

-0.005 

17.00 

5.0to 

0.551 

1.50 

0.009 

7.74 

5.5to 

0.517 

1.70 

0.037 

7.26 

O.Oto 

0.482 

1.89 

0.069 

7.30 



-^thermal “ 

3 

= exp(-t/ri) 

i=l 

U 

ai 

ri(to/h) 

02 

T2{to/h) 

as 

Tsito/h) 

0 

0 

0 

0.631 

9.30 

-0.232 

115.71 

0 


0.5to 

0.631 

10.45 

-0.231 

71.38 

0 


l.Oto 

0.630 

11.74 

-0.310 

61.24 

0.076 

148.85 

1.5to 

0.653 

13.03 

-0.318 

44.19 

0.056 

204.12 

2 . 0 to 

0.649 

13.12 

-0.303 

35.09 

0.042 

264.41 

2.5to 

0.581 

11.73 

- 0.221 

31.21 

0.030 

336.81 

CO 

0 

0 

0.500 

9.61 

- 0.120 

31.34 

0.021 

424.09 

3.5to 

0.461 

7.69 

-0.057 

34.54 

0.014 

530.22 

0 

0 

0.457 

6.27 

-0.026 

40.54 

0.009 

660.50 

4.5to 

0.473 

5.27 

- 0.012 

48.97 

0.006 

820.34 

0 

-Ki 

0 

0.497 

4.58 

-0.006 

57.87 

0.004 

1014.4 

5.5to 

0.522 

4.12 

-0.003 

61.39 

0.003 

1244.2 

6 . 0 to 

0.546 

3.82 

-0.003 

52.30 

0.003 

1511.2 
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